Nonlinear Lattice Relaxation of Photoexcited Diplatinum-Halide Chain Compounds 
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In order to reveal the relaxation mechanism of photogenerated charge-transfer excitations in 
quasi-one-dimensional halogen-bridged diplatinum complexes, we calculate the low-lying adiabatic 
potential energy surfaces of a one-dimensional extended Peierls-Hubbard model. High-energy excita- 
tions above the electron-hole continuum may relax into polarons, while excitons pumped within the 
optical gap are self- localized and then either decay by luminescence or divide into solitons. Neutral 
solitons, charged solitons, and polarons may be simultaneously photogenerated in a diplatinum- 
halide chain, which has never been observed in any conventional platinum-halide chain. Optical 
conductivity is also simulated along the decay paths for experimental verification. 

PACS numbers: 71.45.Lr, 78.20.Bh, 71.35.-y, 78.20.Ci 



I. INTRODUCTION 

Halogen (X)-bridged transition-metal (M) chain com- 
pounds, abbreviated as MX chains, have been attract- 
ing much interest for several decades. Substituting con- 
stituent metals, halogens, ligand molecules, and counter 
ions, we can widely tune their electronic structure 
and systematically study low-dimensional quantum phe- 
nomena in consequence of competing electron-electron 
(e-e) and electron-lattice (e-1) interactions. 1,2 VtX 
chains, including Wolffram's red salt [Pt(C 2 H 7 N)4Cl]Cl2- 
2H2O, 3 exhibit a Peierls-distorted mixed-valent ground 
state, whereas MX chains, which are strongly cor- 
related Mott insulators, have a monovalent regular- 
chain structure. 4 Various ground 5-7 and defect 8-12 
states were theoretically predicted and indeed observed 
experimentally. 13-20 The photoinduced midgap absorp- 
tion of [Pt(C2H 8 N 2 )2Cl] (0104)2 caused a vigorous ar- 
gument on self-trapped nonlinear excitations, 21-24 while 
the gigantic optical nonlinearity of [Ni(CgHi4N2)2Br]Br2 
opened up a new way to optical devices. 25 

The thus-fascinating MX family compounds have been 
gaining renewed interest in recent years due to their bin- 
uclear metal analogs, 26,27 abbreviated as MMX chains, 
which exhibit quantum, 28-30 thermal, 31-33 and pressure- 
induced 34-37 transitions between a wider variety of 
mixed-valent states. 38 ' 39 The direct M(d z 2)-M(d z 2) over- 
lap effectively reduces the on-site Coulomb repulsion and 
therefore makes electrons more itinerant. MMX chains 
are indeed much more conductive than MX chains. 31 
Then charge- and/or spin-carrying local excitations such 
as solitons 40 and polarons 41 are more and more interest- 
ing. Thermally excited spin solitons have already been 
observed in Pt2(C5HnCS2)4l. 42 However, photogener- 
ated MMX defect states have neither been measured nor 
been calculated yet. Effects of metal binucleation on the 
relaxation of photogenerated charge-transfer excitations 
are not only scientifically interesting in themselves but 
must be also the key to optical switching. Thus moti- 
vated, we simulate photoexcitation and nonlinear lattice 
relaxation of diplatinum-halide chains and stimulate fur- 
ther experimental explorations. 



II. CALCULATIONAL PROCEDURE 

We describe MMX chains by the one-dimensional |- 
filled single-band Peierls-Hubbard adiabatic Hamiltonian 

H = - \pMXM - a(/ n +l:- + ^1:+)] ( a n+l, s bn,s 
n,s 

71, S 

-0^2(ln:-n n , s + ln:+m n ,s) + ^ X y^A l n:- 
n,s n 

+£:+)+ 

K,t^n,| + m n ^m n:i ) 

n 

+ J2 ( Vmm 

+ VMXMH n +i >s m„ )S '),(2.1) 



where n n s = oX 



a n>s and m n , s = b} l s b n ^ with a+ s and 
b\ s creating an electron with spin s on the M d z i or- 
bitals in the nth MMX unit. tyiM and tuxu give the in- 
tradimer and intcrdimcr electron hoppings, respectively. 
a and f3 describe the Peierls- and Holstein-type e-1 cou- 
plings, respectively, with Kmx being the metal-halogen 
spring constant. Z„ : _ = v n — u n -\ and l n . + — u n — v n 
with u n and v n being, respectively, the chain-direction 
displacements of the halogen and metal dimer in the 
nth MMX unit from their equilibrium positions. Every 
diplatinum moiety, with its surrounding ligands, is not 
deformed. The notation is further explained in Fig. 1. 
We take tuM as twice tMXM setting tuxu and Kmx both 
equal to unity. Since the parameter sets (a, j3, Kmx) 
and (aa,af3,a 2 KMx) are equivalent to each other with 
an arbitrary constant a, the specific value of Kmx is 
insignificant. 43 The number of MMX units, denoted by 
N, is set equal to 300. 
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FIG. 1: Schematic representation of MMX chains. 
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FIG. 2: (Color online) Adiabatic potential energies as functions of d measured from the ground-state energy, (a): Relaxation 
channels connected to soliton (S)-antisoliton (S) pairs; (b): Relaxation channels connected to polaron (P) pairs. Dotted 
(wavy) arrows suggest possible relaxation paths (luminescence), a — 0.0, (5 = 1.2 and a = 0.3, /? = 0.8 for the P2O5H2- and 
_RCS2-ligand complexes, respectively; Um = 0.5, Vmm = 0.25, and Vmxm = 0.15 in common. 
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FIG. 3: (Color online) Spatial configurations of various photogenerated excitations. We first calculate the local electron 
((i|a^ T an,TlO + ^l a i,4 a n,llO = n "> { l \ b l,-i b n,r\ l ) + { l \ b t,i b n,i\ l ) = m ") and spin ((Z|aJ JjT an,TlO - ('l a ! l ,j. a ™,.|.IO = 2s n , (l\bl A b nA \l)- 
(l\b n |&n,||0 = 2t n ) densities in the Zth excited state and then extract their nonalternating components (n n _i + 2n„ + n n +i = 
-4n„,m„_i + 2m„ + m n +i = —4fh n ; s„-i + 2s n + s n+ i = 4s n ; t n _i + 2t n + t n +\ = At n ). In the above, the average electron 
density is added to the nonalternating components. 
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FIG. 4: (Color online) Adiabatic potential energies as functions of d measured from the ground-state energy, (a): Relaxation 
channels connected to soliton (S)-antisoliton (S) pairs; (b): Relaxation channels connected to polaron (P) pairs. Dotted (wavy) 
arrows suggest possible relaxation paths (luminescence), a = 0.0, f3 = 1.2 and a — 0.3, f3 = 0.8 for the CDW and ACP 
backgrounds, respectively; Um ~ 1-0, Vmm = 0.5, and Vmxm = 0.3 in common. 
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FIG. 5: (Color online) Adiabatic potential energies as func- 
tions of d measured from the ground-state energy in MX 
chains, where we employ the same type of Hamiltonian as 
Eq. (2.1) with a = 0.0, (3 = 0.7, U M = 1.0, and Vmxm = 0.3, 
which is relevant to [Pt^HgNa^XKClO-i^ (X = CI, Br), 
(a): Relaxation channels connected to soliton (S)-antisoliton 
(S) pairs; (b): Relaxation channels connected to polaron (P) 
pairs. Dotted (wavy) arrows suggest possible relaxation paths 
(luminescence). 



The existent diplatinum-halide chain compounds are 
classified into two groups: At[Pt 2 (P 2 05H 2 )4A]-nH 2 
(X = CI, Br, I; A = Li, Cs, • • -) 26 ' 44 and Pt 2 (i?CS 2 ) 4 I 
(R = C„H 2n+ i). 27 ' 45 The former structurally re- 
sembles the MX conventional and exhibits a 
ground state with halogen-sublattice dimerization: 

-X-- ■ ■ Pt 2+ Pt 2+ - • ■ X- - Pt 3+ Pt 3+ - X~- ■ ; Which 

is referred to as the charge-density-wave (CDW) state, 
where the intrasite e-1 coupling fj is dominant. The latter 
possesses a distinct ground state with metal-sublattice 
dimerization, - I"- • • Pt 2+ Pt 3+ -r-Pt 3+ Pt 2+ - • -I"- • •, 
which is referred to as the alternate charge-polarization 
(ACP) state, where the intersite e-1 coupling a is 
significant. Platinum-halide chains exhibit intermediate 
e-e interactions (Vmxm *C Um % £mxm) 2,9 and platinum 
binucleation should reduce the on-site repulsion. Thus 
the Pt 2 X Coulomb parameters are, unless otherwise 
noted, set for Um — 0.5, Vmm — 0.25, and Vmxm = 0.15. 
The e-1 coupling constants are taken in two ways as 
a = 0.0, f3 = 1.2 and a — 0.3, (3 = 0.8, which are relevant 
to the P 2 0sH 2 - and i?CS 2 -ligand complexes and indeed 
give the CDW and ACP ground states, respectively, 
under the above Coulomb parametrization. 46 

Photogenerated charge-transfer excitations, spreading 
over the chain at first, are self- localized into excitons, 
solitons, and polarons. The whole relaxation scenario is 
describable with a trial wave function 9,41 ' 47 



»n:± 



-l) n /o 



1 + 511 tanh 



\n±Sn\ -d/2 



. (2-2) 



where lo is set equal to the halogen-ion displacement in 
the uniformly distorted CDW ground state, while the rest 
of the variational parameters are determined so as to min- 
imize the energy of the lowest-lying excited state. Once 
the CDW state is photoexcited into the Frank-Condon 
state, which still sits at 51 = 0, the uniform bond alter- 
nation begins to be locally deformed. Increasing 51 with d 
fixed to zero represents the self-localization of a charge- 
transfer exciton. The self-trapped exciton (STE) may 



further divide into a pair of local defects with increasing 
d. There is also a possibility of a higher-energy pumped 
electron-hole (e-h) pair directly splitting into distant de- 
fects in a pair, with simultaneously increasing 51 and d. 
£ corresponds to the extent of a local defect, whereas 5n 
allows the neighboring metal sites to behave unequally 
around the defect center. <r„ takes ± according to relax- 
ation channels. It is simply set equal to unity for solitonic 
defects, while it is defined as a n — sgn(|n ± 5n\ — d/2) 
for polaronic defects. In both cases, the key variable d 
indicates the interdefect distance. We optimize 51, £, and 
5n at each d. 

With given l n -.±, we solve the Hartree-Fock (HF) 
Hamiltonian 7Yhf and obtain the low-lying states |£)hf 
of energy Ef^ (/ = 0, 1, 2, • • •). Then we consider refining 
the description of excited states, which may be generally 
given as 



HF' 



(2.3) 



where ep is the Fermi energy and c\ creates an elec- 
tron with spin s in the Ath HF eigenstate of energy e\ >s . 
The HF scheme describes any excited state as a single 
Slater determinant, that is, f(/j,,i/,s;l) — 5^ s .i, taking 
no account of the residual interaction H — Hhf = V. 
Full diagonalization of TL on the basis of \l) must be 
the best way to deal with the excitonic effect. How- 
ever, even at = 120, for instance, such a calculation 
costs eight gigabyte memory and a hundred hours, which 
means spending more than a year in optimizing the wave 
function (2.2) at every fixed d. In an attempt to take 
the excitonic effect into calculation more efficiently, we 
may consider neglecting all the off-diagonal elements of 
V and correcting the energy scheme perturbationally. 9 
Then the Ith excited-state energy is expressed as Ei = 
Eq + e VtS - e^s + hf^MOhfi wnere E o = E^ F . Such 
a perturbational treatment of the excitonic effect is not 
only practical but also fairly quantitative under the not- 
so-strong electronic correlation of our present interest. 
We have indeed confirmed for short chains that the thus- 
obtained energies E\ well approximate those of excited 
states of the configuration-interaction (CI) type (see Ap- 
pendix A). This is not the case with strongly correlated 
nickel complexes. 48 The following variational calculation 
is carried out so as to minimize E\ with respect to 51, 
£, and 5n at each d. When we take particular interest 
in the ith energy surface, Ei may be minimized instead. 
However, the whole energy scheme is not sensitive to the 
variational target and remains unchanged visually. 



III. ADIABATIC POTENTIAL ENERGY 
SURFACES 

Figure 2 presents the thus-calculated energy 
surfaces, the left two of which are relevant to 
^ 4 [Pt 2 (P 2 5 H 2 )4A]-nH 2 0, while the right two of 
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which arc to Pt 2 (i?CS 2 )4l. The electron-hole contin- 
uum of Pt 2 (i?CS 2 )4l lies lower in energy than that 
of A 4 [Pt 2 (P 2 5 H 2 ) 4 A]-nH 2 0. In fact Pt 2 (CH 3 CS 2 )4l 
exhibits much smaller optical gap than conventional MX 
chain compounds and even shows metallic conduction 
at room temperature. 31 The relaxation channels of the 
P 2 C>5H 2 - and i?CS 2 -ligand complexes are qualitatively 
the same. Fully trapped excitons [Fig. 3(a)], unless 
decay by luminescence, are dissociated into soliton 
(S)-antisoliton (S) pairs [Figs. 3(b) and 3(c)], whereas 
there is no possibility of their relaxing into polaron 
(P) pairs. P~^ — P + ^ pairs [Fig. 3(d)] may be created 
from higher-energy excited states. No energy barrier 
between an STE and any S — S pair, which is not the 
case with PtX chains, 9 ' 47 should result in very short 
decay time of STEs and long life time of S — S pairs. The 
soliton and antisoliton in any S — S pair have localized 
wave functions and their overlap rapidly decreases with 
increasing d. Therefore, instantaneous charge transport 
between far distant S and S is hardly probable and 
tunneling between the energy surfaces E\ and Eq is 
restricted to the region of moderately small d. 

Under the relatively weak correlation relevant to Pt 2 A 
chains, charged soliton pairs S~ — S + are lower in en- 
ergy than neutral soliton pairs S ^ — S°^. 49 Figures 3(b) 
and 3(c) show that solitons and antisolitons on the en- 
ergy surface E\ carry net spins, whereas those on the 
energy surface Eq convey net charges. As the Coulomb 
repulsion grows, the polaronic channel qualitatively re- 
mains unchanged, while the solitonic one significantly 
varies. Figure 4 indeed shows that S~ — S + pairs are 
possibly higher in energy than STEs as well as — S°^ 
pairs. Increasing on-site Coulomb repulsion reduces the 
Peierls gap, especially on the CDW background. Fig- 
ure 4 is reminiscent of the situation in VtX chains, 9 ' 47 
which is reproduced in Fig. 5. Photoinduccd midgap 
absorption of [Pt(C 2 H 8 N 2 ) 2 X] (C10 4 ) 2 (X = CI, Br) 18 ' 20 
was actually attributed to neutral solitons. 4,47 In com- 
parison with PtCl and PtBr chains, PtI chains have a 
larger supertransfcr energy tuxu (Ref. 1) and thus 
exhibit effectively suppressed e-e interactions. Conse- 
quently there appear charged solitons instead of neu- 
tral ones in photoexcitcd [Pt(C 2 H 8 N 2 ) 2 I] (C10 4 ) 2 - 17 Pt-X" 
compounds generally exhibit further absorption bands 
due to polarons within the gap, 4 ' 50 provided the exci- 
tation energy is higher than the charge-transfer gap A 
However, there is no signal of neutral and charged soli- 
tons being simultaneously photogenerated in PtA chains, 
probably because either S 0J -— S° T or S~— S + pairs are nec- 
essarily higher in energy than STEs. 47 Figure 2 promises 
dramatic observations: Neutral solitons, charged soli- 
tons, and polarons may be simultaneously detected in 
photoexcitcd diplatinum-halide chain compounds. The 
decay time of luminescent STEs, r, is a few hundred 
picoseconds in [Pt(C 2 H 8 N 2 ) 2 Cl] (C10 4 ) 2 , 51 ' 52 which is 
much shorter than the radiative lifetime and therefore 
comes from the dissociation to S ^ — S° T pairs. In 



At[Pt 2 (P 2 5 H 2 ) 4 A]-nH 2 0, the S oi - S° T path is likely 
open to an STE without any potential barrier in between 
and therefore STEs must decay still faster possibly with 
r < 100 ps. 



IV. PHOTOINDUCED ABSORPTION 
SPECTRA 

In order to encourage time-resolved optical measure- 
ments on Pt 2 A chains, we calculate photoinduced ab- 
sorption spectra as functions of the interdefect distance 
d. The real part of the optical conductivity on the po- 
tential energy surface E\ is represented as 

ff H = £| WI0| 2 W -Ei- M, (4-1) 

where the current operator J = ^ n g j„, s is given by 

jn,s — -J^CMXM \pMXM - a(l n +l:- + ln:+)] (<4i+l, s &n,s 

-bl s a n+1 , s ) + -c M MtMM(bi ),(4-2) 

with cmm and Cmxm being the average M-M and M- 
X-M distances, respectively, and set for Cmxm — 2cmm- 
When we take |Z)hf f° r 10 with the perturbationally cor- 
rected energy scheme Ei, we obtain \l'} — c' v , s /<Vv |Z) HF 

and E v = Ef v + t v , , s , - e^.y + (l'\V\l'). Figures 
6(a)— 6(d) present the optical conductivity spectra along 
the relaxation paths labeled 1—4 in Fig. 2, where (a) — (c) 
are characteristic of solitonic excitations via STEs, while 
(d) of a polaronic one from the e-h continuum. The intra- 
gap absorption bands visualize various nonlinear lattice 
relaxation paths distinguishably, which are scaled up and 
interpreted in Fig. 7. 

In Fig. 7(a), an STE splits into far distant S ^ and S°\ 
while in Fig. 7(b), an STE relaxes into far distant S~ 
and S + via quantum tunneling. In PtX chains, S~— S + 
pairs, instead of S oi -S T pairs, lie in the same potential 
energy surface as STEs', but they can not be reached 
due to their high energy. Here in Pt 2 X chains, STEs 
can nonradiatively decay into both neutral- and charged- 
soliton pairs and therefore the early decrease of their 
luminescence intensity should look double-exponential. 
The optical conductivity spectra are informative about 
the soliton charge. With increasing d, the exciton-to- 
band infrared absorption peak moves upward, whereas 
the intra-exciton midgap one downward. After passing 
each other, the former merges into the soliton-to-band 
absorption, while the latter fades out with S and S going 
away from and less overlapping with each other. Such a 
scenario is common to both S ^— S°^ and S~— S + paths, 
but their long-time spectra look different. The neutral- 
soliton-to-band absorption spectrum consists of two close 
peaks and is thus broader than the charged-soliton-to- 
band one. For Pt 2 (i?CS 2 )4l, it is remarkably humped 
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FIG. 6: (Color online) Photoinduced absorption spectra along possible relaxation paths: (a) STE — > S ^ — S° T ; (b) STE — > 
S" - S + ; (c) STE -> S - S -<■ CDW; (d) e - h -> P~ T - P +i . 



due to the two intragap soliton levels lying relatively 
far apart from each other. On the ACP background, 
the intragap soliton levels are more sensitive to the e- 
e interactions. 53 Pt2(i?CS2)4l exhibits another surviving 
absorption in the infrared region, which is attributable to 
the electron transfer between S and S . This is because of 
larger £ in Pt 2 (i?CS 2 )4l (see Fig. 3). Solitons are more 
delocalized on the ACP background than on the CDW 
background. Figure 7(c) describes a nonradiative gemi- 
nate recombination through the solitonic state, where the 
intra-exciton absorption peak turns upward and merges 
with the background intcrband absorption. 

Figure 7(d) depicts polaronic excitations, which are ef- 
ficiently generated through high-energy pumping. There 
are four potential energy surfaces leading to P~ T — P + ^ 
pairs. Although STEs lie in one of them, they have no 
chance of relaxing into polarons. When charge-transfer 
excitations split into polaron pairs, the intragap absorp- 
tion spectra at small d, that is, the early observations, 
significantly vary with their relaxation paths, but they all 
converge to well-separated two bands. The intra-polaron 
absorption is higher in energy and stronger in intensity 
than the polaron-to-band one. 



V. SUMMARY 

We have calculated the nonlinear lattice relaxation 
paths of photogenerated charge-transfer excitations in 
diplatinum-halide chain compounds and "observed" their 
optical conductivity spectra. There is a possibility of 
neutral solitons, charged solitons, and polarons coexist- 
ing in photoexcited Pt^X chains especially of the i?CS2- 
ligand type, which is never the case with conventional 
PtX chains. The optical-conductivity spectra due to po- 



laron pairs generally comprise well-separated two bands, 
while those due to soliton pairs mainly consist of a single 
band but its structure is rather varied with the ground 
state and the relaxation path. The S~-to-band absorp- 
tion is truly single-peaked, while the S°-to-band absorp- 
tion gives a broader band with a shoulder in the lower- 
energy side. In Pt2(-ffCS2)4l, the S-to-S absorption is 
long surviving and detectable in the infrared region in 
addition. 

Indeed there are many similarities between the Pt 2 X 
and PtA relaxation mechanisms, but their solitonic chan- 
nels qualitatively differ. In Pt^X chains, both neutral- 
and charged- soliton paths are open to photogenerated 
charge-transfer excitons and soliton pairs stably sur- 
vive longer life time. The present study is motivated 
in part by the pioneering calculations on PtX chains, 9 
where no observable to be measured was presented, how- 
ever. Nowadays, the time-resolved optical spectroscopy 
technique has been much more refined than before and 
the full spectrum can be analyzed with femtosecond 
precision. 54 Such circumstances have stimulated us to 
have the idea of calculating Fig. 6, spending a week 
on each run along the relaxation paths. Even the pho- 
toinduced midgap absorption of a prototypical platinum- 
halide chain compound [Pt(C 2 H 8 N 2 )2Cl](C104)2 14 ' 55 is 
not yet completely solved in spite of the enthusias- 
tic argument. 15 ' 18,22 ' 23,47 Diplatinum-halide chain com- 
pounds comprise two groups with distinct mixed-valent 
ground states and their electronic states are much more 
tunable. Comparative measurements, including photoin- 
duced absorption, luminescence, and electron spin res- 
onance, on the P2O5H2- and i?CS2-ligand MMX com- 
plexes must reveal their intrinsic nonlinear lattice relax- 
ation mechanism and may even give a key to unsettled 
issues on MX complexes. 
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FIG. 7: (Color online) The low-energy part of Fig. 6 is scaled up, together with the intragap level scheme and occupancy, where 
round arrows denote optically allowed excitations, while dotted arrows suggest tunneling between different energy surfaces. 
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with Figs. 2(a) and 4(a) in Fig. 8. 
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APPENDIX A: COMPARISON BETWEEN THE 
HF AND CI SCHEMES 

Adiabatic potential energies of short chains (N = 50) 
are calculated by the CI scheme (x) as well as the per- 



7 



CDW background ACP background 




10 20 30 10 20 30 



FIG. 8: (Color online) Relaxation paths leading to soliton (S)- 
antisoliton (S) pairs, a = 0.0, f3 = 1.2 and a = 0.3, f3 — 0.8 for 
the CDW and ACP backgrounds, respectively, (a) Um = 0.5, 
Vmm = 0.25, Vmxm = 0.15; (b) U M = 1.0, Vmm = 0.5, 
Vmxm =0.3. 
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